Brine-driven destruction of clay minerals in Gale crater, Mars.

Mars' sedimentary rock record preserves information on geological (and potential astrobiological) processes that occurred on the planet billions of years ago. The Curiosity rover is exploring the lower reaches of Mount Sharp, in Gale crater on Mars. A traverse from Vera Rubin ridge to Glen Torridon has allowed Curiosity to examine a lateral transect of rock strata laid down in a martian lake ~3.5 billion years ago. We report spatial differences in the mineralogy of time-equivalent sedimentary rocks <400 meters apart. These differences indicate localized infiltration of silica-poor brines, generated during deposition of overlying magnesium sulfate-bearing strata. We propose that destabilization of silicate minerals driven by silica-poor brines (rarely observed on Earth) was widespread on ancient Mars, because sulfate deposits are globally distributed.

Evidence for Multiple Diagenetic Episodes in Ancient Fluvial-Lacustrine Sedimentary Rocks in Gale Crater, Mars.

The Curiosity rover's exploration of rocks and soils in Gale crater has provided diverse geochemical and mineralogical data sets, underscoring the complex geological history of the region. We report the crystalline, clay mineral, and amorphous phase distributions of four Gale crater rocks from an 80-m stratigraphic interval. The mineralogy of the four samples is strongly influenced by aqueous alteration processes, including variations in water chemistries, redox, pH, and temperature. Localized hydrothermal events are evidenced by gray hematite and maturation of amorphous SiO2 to opal-CT. Low-temperature diagenetic events are associated with fluctuating lake levels, evaporative events, and groundwater infiltration. Among all mudstones analyzed in Gale crater, the diversity in diagenetic processes is primarily captured by the mineralogy and X-ray amorphous chemistry of the drilled rocks. Variations indicate a transition from magnetite to hematite and an increase in matrix-associated sulfates suggesting intensifying influence from oxic, diagenetic fluids upsection. Furthermore, diagenetic fluid pathways are shown to be strongly affected by unconformities and sedimentary transitions, as evidenced by the intensity of alteration inferred from the mineralogy of sediments sampled adjacent to stratigraphic contacts.

Chemistry, mineralogy, and grain properties at Namib and High dunes, Bagnold dune field, Gale crater, Mars: A synthesis of Curiosity rover observations.

The Mars Science Laboratory Curiosity rover performed coordinated measurements to examine the textures and compositions of aeolian sands in the active Bagnold dune field. The Bagnold sands are rounded to subrounded, very fine to medium sized (~45-500 µm) with ≥6 distinct grain colors. In contrast to sands examined by Curiosity in a dust-covered, inactive bedform called Rocknest and soils at other landing sites, Bagnold sands are darker, less red, better sorted, have fewer silt-sized or smaller grains, and show no evidence for cohesion. Nevertheless, Bagnold mineralogy and Rocknest mineralogy are similar with plagioclase, olivine, and pyroxenes in similar proportions comprising >90% of crystalline phases, along with a substantial amorphous component (35% ± 15%). Yet Bagnold and Rocknest bulk chemistry differ. Bagnold sands are Si enriched relative to other soils at Gale crater, and H2O, S, and Cl are lower relative to all previously measured Martian soils and most Gale crater rocks. Mg, Ni, Fe, and Mn are enriched in the coarse-sieved fraction of Bagnold sands, corroborated by visible/near-infrared spectra that suggest enrichment of olivine. Collectively, patterns in major element chemistry and volatile release data indicate two distinctive volatile reservoirs in Martian soils: (1) amorphous components in the sand-sized fraction (represented by Bagnold) that are Si-enriched, hydroxylated alteration products and/or H2O- or OH-bearing impact or volcanic glasses and (2) amorphous components in the fine fraction (<40 µm; represented by Rocknest and other bright soils) that are Fe, S, and Cl enriched with low Si and adsorbed and structural H2O.

The first X-ray diffraction measurements on Mars.

The Mars Science Laboratory landed in Gale crater on Mars in August 2012, and the Curiosity rover then began field studies on its drive toward Mount Sharp, a central peak made of ancient sediments. CheMin is one of ten instruments on or inside the rover, all designed to provide detailed information on the rocks, soils and atmosphere in this region. CheMin is a miniaturized X-ray diffraction/X-ray fluorescence (XRD/XRF) instrument that uses transmission geometry with an energy-discriminating CCD detector. CheMin uses onboard standards for XRD and XRF calibration, and beryl:quartz mixtures constitute the primary XRD standards. Four samples have been analysed by CheMin, namely a soil sample, two samples drilled from mudstones and a sample drilled from a sandstone. Rietveld and full-pattern analysis of the XRD data reveal a complex mineralogy, with contributions from parent igneous rocks, amorphous components and several minerals relating to aqueous alteration. In particular, the mudstone samples all contain one or more phyllosilicates consistent with alteration in liquid water. In addition to quantitative mineralogy, Rietveld refinements also provide unit-cell parameters for the major phases, which can be used to infer the chemical compositions of individual minerals and, by difference, the composition of the amorphous component.

Mineralogy of a mudstone at Yellowknife Bay, Gale crater, Mars.

Sedimentary rocks at Yellowknife Bay (Gale crater) on Mars include mudstone sampled by the Curiosity rover. The samples, John Klein and Cumberland, contain detrital basaltic minerals, calcium sulfates, iron oxide or hydroxides, iron sulfides, amorphous material, and trioctahedral smectites. The John Klein smectite has basal spacing of ~10 angstroms, indicating little interlayer hydration. The Cumberland smectite has basal spacing at both ~13.2 and ~10 angstroms. The larger spacing suggests a partially chloritized interlayer or interlayer magnesium or calcium facilitating H2O retention. Basaltic minerals in the mudstone are similar to those in nearby eolian deposits. However, the mudstone has far less Fe-forsterite, possibly lost with formation of smectite plus magnetite. Late Noachian/Early Hesperian or younger age indicates that clay mineral formation on Mars extended beyond Noachian time.

X-ray diffraction results from Mars Science Laboratory: mineralogy of Rocknest at Gale crater.

The Mars Science Laboratory rover Curiosity scooped samples of soil from the Rocknest aeolian bedform in Gale crater. Analysis of the soil with the Chemistry and Mineralogy (CheMin) x-ray diffraction (XRD) instrument revealed plagioclase (~An57), forsteritic olivine (~Fo62), augite, and pigeonite, with minor K-feldspar, magnetite, quartz, anhydrite, hematite, and ilmenite. The minor phases are present at, or near, detection limits. The soil also contains 27 ± 14 weight percent x-ray amorphous material, likely containing multiple Fe(3+)- and volatile-bearing phases, including possibly a substance resembling hisingerite. The crystalline component is similar to the normative mineralogy of certain basaltic rocks from Gusev crater on Mars and of martian basaltic meteorites. The amorphous component is similar to that found on Earth in places such as soils on the Mauna Kea volcano, Hawaii.

Curiosity at Gale crater, Mars: characterization and analysis of the Rocknest sand shadow.

The Rocknest aeolian deposit is similar to aeolian features analyzed by the Mars Exploration Rovers (MERs) Spirit and Opportunity. The fraction of sand <150 micrometers in size contains ~55% crystalline material consistent with a basaltic heritage and ~45% x-ray amorphous material. The amorphous component of Rocknest is iron-rich and silicon-poor and is the host of the volatiles (water, oxygen, sulfur dioxide, carbon dioxide, and chlorine) detected by the Sample Analysis at Mars instrument and of the fine-grained nanophase oxide component first described from basaltic soils analyzed by MERs. The similarity between soils and aeolian materials analyzed at Gusev Crater, Meridiani Planum, and Gale Crater implies locally sourced, globally similar basaltic materials or globally and regionally sourced basaltic components deposited locally at all three locations.

Bond paths and van der Waals interactions in orpiment, As2S3.

The calculated electron density distribution for orpiment, As(2)S(3), reveals that As-S, S-S, and As-As bond paths are associated with the experimental interlayer directed bonded interactions detected in a combined infrared and Raman study. The successful modeling of the infrared- and Raman-determined interlayer bonded interactions together with bond paths and the structuralization of a variety of inorganic molecules in terms of "key-lock" bond path mainstays support the argument that van der Waals forces in inorganic molecular crystals are directional.

Role of directed van der Waals bonded interactions in the determination of the structures of molecular arsenate solids.

Bond paths, local energy density properties, and Laplacian, L(r) = -wedge(2)rho(r), composite isosurfaces of the electron density distributions were calculated for the intramolecular and intermolecular bonded interactions for molecular solids of As(2)O(3) and AsO(2) composition, an As(2)O(5) crystal, a number of arsenate molecules, and the arsenic metalloid, arsenolamprite. The directed intermolecular van der Waals As-O, O-O, and As-As bonded interactions are believed to serve as mainstays between the individual molecules in each of the molecular solids. As-O bond paths between the bonded atoms connect Lewis base charge concentrations and Lewis acid charge depletion domains, whereas the O-O and As-As paths connect Lewis base pair and Lewis acid pair domains, respectively, giving rise to sets of intermolecular directed bond paths. The alignment of the directed bond paths results in the periodic structures adopted by the arsenates. The arrangements of the As atoms in the claudetite polymorphs of As(2)O(3) and the As atoms in arsenolamprite are similar. Like the As(2)O(3) polymorphs, arsenolamprite is a molecular solid connected by relatively weak As-As intermolecular directed van der Waals bond paths between the layers of stronger As-As intramolecular bonded interactions. The bond critical point and local energy density properties of the intermolecular As-As bonded interactions in arsenolamprite are comparable with the As-As interactions in claudetite I. As such, the structure of claudetite I can be viewed as a stuffed derivative of the arsenolamprite structure with O atoms between pairs of As atoms comprising the layers of the structure. The cubic structure adopted by the arsenolite polymorph can be understood in terms of sets of directed acid-base As-O and base-base O-O pair domains and bond paths that radiate from the tetrahedral faces of its constituent molecules, serving as face-to-face key-lock mainstays in forming a periodic tetrahedral array of molecules rather than one based on some variant of close packing. The relatively dense structure and the corrugation of the layers in claudetite I can also be understood in terms of directed van der Waals As-O bonded interactions. Our study not only provides a new basis for understanding the crystal chemistry and the structures of the arsenates, but it also calls for a reappraisal of the concept of van der Waals bonded interactions, how the structures of molecular solids are viewed, and how the molecules in these solids are bonded in a periodic structure.

Raman spectroscopic and visible absorption investigation of LiCrSi2O6 pyroxene under pressure.

The first observation of the vibrational spectrum of the synthetic pyroxene Li-kosmochlor (LiCrSi2O6) is reported herein. The Raman and visible spectra are reported as a function of pressure. Though the pyroxene retains its P21/c symmetry, changes in the Raman spectra are observed between 6.8 and 7.7 GPa, possibly due to the formation of an additional bond between Li and O3 or some other transition that retains the mineral's P21/c space group. Splitting of the peak appearing at approximately 700 cm(-1), used to characterize the P21/c phase in other studies, is not observed. Comparison is made with the Raman spectra of LiAlSi2O6 and LiFeSi2O6 in the P21/c phase and the visible spectra of NaCrSi2O6 at high pressures.

Experimental bond critical point and local energy density properties determined for Mn-O, Fe-O, and Co-O bonded interactions for tephroite, Mn2SiO4, fayalite, Fe2SiO4, and Co2SiO4 olivine and selected organic metal complexes: comparison with properties calculated for non-transition and transition metal M-O bonded interactions for silicates and oxides.

Bond critical point (bcp) and local energy density properties for the electron density (ED) distributions, calculated with first-principle quantum mechanical methods for divalent transition metal Mn-, Co-, and Fe-containing silicates and oxides are compared with experimental model ED properties for tephroite, Mn 2SiO 4, fayalite, Fe 2SiO 4, and Co 2SiO 4 olivine, each determined with high-energy synchrotron single-crystal X-ray diffraction data. Trends between the experimental bond lengths, R(M-O), (M = Mn, Fe, Co), and the calculated bcp properties are comparable with those observed for non-transition M-O bonded interactions. The bcp properties, local total energy density, H( r c), and bond length trends determined for the Mn-O, Co-O, and Fe-O interactions are also comparable. A comparison is also made with model experimental bcp properties determined for several Mn-O, Fe-O, and Co-O bonded interactions for selected organometallic complexes and several oxides. Despite the complexities of the structures of the organometallic complexes, the agreement between the calculated and model experimental bcp properties is fair to good in several cases. The G( r c)/rho( r c) versus R(M-O) trends established for non-transition metal M-O bonded interactions hold for the transition metal M-O bonded interactions with G( r c)/rho( r c) increasing in value as H( r c) becomes progressively more negative in value, indicating an increasing shared character of the interaction as G( r c)/rho( r c) increases in value. As observed for the non-transition metal M-O bonded interactions, the Laplacian, nabla (2)rho( r c), increases in value as rho( r c) increases and as H( r c) decreases and becomes progressive more negative in value. The Mn-O, Fe-O, and Co-O bonded interactions are indicated to be of intermediate character with a substantial component of closed-shell character compared with Fe-S and Ni-S bonded interactions, which show greater shared character based on the | V( r c)|/ G( r c) bond character indicator. The atomic charges conferred on the transition metal atoms for the three olivines decrease with increasing atomic number from Mn to Fe to Co as the average M-O bond lengths decrease from 2.219 to 2.168 to 2.128 A, respectively.

Shared and closed-shell O-O interactions in silicates.

Bond paths of maximum electron density spanning O-O edges shared between equivalent or quasiequivalent MOn (n > 4) coordination polyhedra are not uncommon electron density features displayed by silicates. On the basis of the positive values for the local electronic energy density, H(rc), at the bond critical points, rc, they qualify as weak "closed-shell" interactions. As observed for M-O bonded interactions (M = first and second row metal atoms), the electron density, rho(rc), and the Laplacian of the electron density increase in a regular way as the separation between the O atoms, R(O-O), decreases. A simple model, based on R(O-O) and the distances of the Si atoms from the midpoint between adjacent pairs of O atoms, partitions the O-O bond paths in the high-pressure silica polymorph coesite into two largely disjoint domains, one with and one without bond paths. The occurrence of O-O bond paths shared in common between equivalent coordination polyhedra suggests that they may be grounded in some cases on factors other than bonded interactions, particularly since they are often displayed by inert procrystal representations of the electron density. In these cases, it can be argued that the accumulation of the electron density along the paths has its origin, at least in part, in the superposition of the peripheral electron density distributions of the metal M atoms occupying the edge-sharing polyhedra. On the other hand, the accumulation of electron density along the paths may stabilize a structure by shielding the adjacent M atoms in the edge-sharing polyhedra. For closed-shell Li-O, Na-O, and Mg-O interactions, H(rc) is positive and increases as the value of rho(rc) increases, unlike the "shared" Be-O, B-O, C-O, Al-O, Si-O, P-O, and S-O interactions, where H(rc) is negative and decreases as rho(rc) increases. The H(rc) values for the weak closed-shell O-O interactions also increase as rho(rc) increases, as observed for the closed-shell M-O interactions. On the basis of the bond critical point properties and the negative H(rc) value, the O-O interaction comprising the O2 molecule in silica III qualifies as a shared interaction.

Theoretical electron density distributions for Fe- and Cu-sulfide earth materials: a connection between bond length, bond critical point properties, local energy densities, and bonded interactions.

Bond critical point and local energy density properties together with net atomic charges were calculated for theoretical electron density distributions, rho(r), generated for a variety of Fe and Cu metal-sulfide materials with high- and low-spin Fe atoms in octahedral coordination and high-spin Fe atoms in tetrahedral coordination. The electron density, rho(rc), the Laplacian, triangle down2rho(rc), the local kinetic energy, G(rc), and the oxidation state of Fe increase as the local potential energy density, V(rc), the Fe-S bond lengths, and the coordination numbers of the Fe atoms decrease. The properties of the bonded interactions for the octahedrally coordinated low-spin Fe atoms for pyrite and marcasite are distinct from those for high-spin Fe atoms for troilite, smythite, and greigite. The Fe-S bond lengths are shorter and the values of rho(rc) and triangle down2rho(rc) are larger for pyrite and marcasite, indicating that the accumulation and local concentration of rho(r) in the internuclear region are greater than those involving the longer, high-spin Fe-S bonded interactions. The net atomic charges and the bonded radii calculated for the Fe and S atoms in pyrite and marcasite are also smaller than those for sulfides with high-spin octahedrally coordinated Fe atoms. Collectively, the Fe-S interactions are indicated to be intermediate in character with the low-spin Fe-S interactions having greater shared character than the high-spin interactions. The bond lengths observed for chalcopyrite together with the calculated bond critical point properties are consistent with the formula Cu+Fe3+S2. The bond length is shorter and the rho(rc) value is larger for the FeS4 tetrahedron displayed by metastable greigite than those displayed by chalcopyrite and cubanite, consistent with a proposal that the Fe atom in greigite is tetravalent. S-S bond paths exist between each of the surface S atoms of adjacent slabs of FeS6 octahedra comprising the layer sulfide smythite, suggesting that the neutral Fe3S4 slabs are linked together and stabilized by the pathways of electron density comprising S-S bonded interactions. Such interactions not only exist between the S atoms for adjacent S8 rings in native sulfur, but their bond critical point properties are similar to those displayed by the metal sulfides.

Classification of metal-oxide bonded interactions based on local potential- and kinetic-energy densities.

A classification of the hydrogen fluoride H-F-bonded interactions comprising a large number of molecules has been proposed by Espinosa et al. [J. Chem. Phys. 117, 5529 (2002)] based on the ratio /Vr(c)/ / Gr(c) where /Vr(c)/ is the magnitude of the local potential-energy density and Gr(c) is the local kinetic-energy density, each evaluated at a bond critical point r(c). A calculation of the ratio for the M-O bonded interactions comprising a relatively large number of oxide molecules and earth materials, together with the constraints imposed by the values of inverted Delta2rho r(c) and the local electronic energy density, Hr(c) = Gr(c) + Vr(c), in the H-F study, yielded practically the same classification for the oxides. This is true despite the different trends that hold between the bond critical point and local energy density properties with the bond lengths displayed by the H-F and M-O bonded interactions. On the basis of the ratio, Li-O, Na-O, and Mg-O bonded interactions classify as closed-shell ionic bonds, Be-O, Al-O, Si-O, B-O, and P-O interactions classify as bonds of intermediate character with the covalent character increasing from Be-O to P-O. N-O interactions classify as shared covalent bonds. C-O and S-O bonded interactions classify as both intermediate and covalent bonded interactions. The C-O double- and triple-bonded interactions classify as intermediate-bonded interactions, each with a substantial component of covalent character and the C-O single-bonded interaction classifies as a covalent bond whereas their local electronic energy density values indicate that they are each covalent bonded interactions. The ratios for the Be-O, Al-O, and Si-O bonded interactions indicate that they have a substantial component of ionic character despite their classification as bonds of intermediate character. The trend between the ratio and the character of the bonded interactions is consistent with trends expected from electronegativity considerations. The ratio increases as the net charges and the coordination numbers for the atoms for several Ni-sulfides decrease. On the contrary, the ratio for the Si-O bonded interactions for the orthosilicate, forsterite, Mg2SiO4, and the high-pressure silica polymorph, stishovite, decreases as the observed net atomic charges and the coordination numbers of Si and O increase in value. The ratio for the Ni-Ni bonded interactions for the Ni-sulfides and bulk Ni metal indicate that the interactions are intermediate in character with a substantial component of ionic character.

Comparison of the electron localization function and deformation electron density maps for selected earth materials.

The electron localization function (ELF) and experimental and theoretical deformation electron density maps are compared for several earth materials and one representative molecule. The number and arrangement of the localized one-electron probability density domains generated in a mapping of the ELF correspond to the number and arrangement of the localized electron density domains generated in a mapping of the deformation electron density distribution, a correspondence that suggests that the two fields are homeomorphically related. As a homeomorphic relationship has been established previously between the Laplacian of the electron density distribution and the ELF, the relationship suggests that the deformation electron density distribution is also homeomorphically related to the Laplacian of the distribution.

Electron density distributions calculated for the nickel sulfides millerite, vaesite, and heazlewoodite and nickel metal: a case for the importance of ni-ni bond paths for electron transport.

Bond paths and the bond critical point properties (the electron density (rho) and the Hessian of rho at the bond critical points (bcp's)) have been calculated for the bonded interactions comprising the nickel sulfide minerals millerite, NiS, vaesite, NiS(2), and heazlewoodite, Ni(3)S(2), and Ni metal. The experimental Ni-S bond lengths decrease linearly as the magnitudes of the properties each increases in value. Bond paths exist between the Ni atoms in heazlewoodite and millerite for the Ni-Ni separations that match the shortest separation in Ni metal, an indicator that the Ni atoms are bonded. The bcp properties of the bonded interactions in Ni metal are virtually the same as those in heazlewoodite and millerite. Ni-Ni bond paths are absent in vaesite where the Ni-Ni separations are 60% greater than those in Ni metal. The bcp properties for the Ni-Ni bonded interactions scatter along protractions of the Ni-S bond length-bcp property trends, suggesting that the two bonded interactions have similar characteristics. Ni-Ni bond paths radiate throughout Ni metal and the metallic heazlewoodite structures as continuous networks whereas the Ni-Ni paths in millerite, a p,d-metal displaying ionic and covalent features, are restricted to isolated Ni(3) rings. Electron transport in Ni metal and heazlewoodite is pictured as occurring along the bond paths, which behave as networks of atomic size wires that radiate in a contiguous circuit throughout the two structures. Unlike heazlewoodite, the electron transport in millerite is pictured as involving a cooperative hopping of the d-orbital electrons from the Ni(3) rings comprising Ni(3)S(9) clusters to Ni(3) rings in adjacent clusters via the p-orbitals on the interconnecting S atoms. Vaesite, an insulator at low temperatures and a doped semiconductor at higher temperatures, lacks Ni-Ni bond paths. The net charges conferred on the Ni and S atoms are about a quarter of their nominal charges for the atoms in millerite and vaesite with the net charge on Ni increasing with increasing Ni-S bond length. Reduced net charges are observed on the Ni atoms in heazlewoodite and are related to its Ni-Ni metal bonded interactions and to the greater covalent character of its bonds. Local energy density and bond critical point properties of the electron density distributions indicate that the Ni-S and Ni-Ni bonded interactions are intermediate in character between ionic and covalent.

Systematic generation of all nonequivalent closest-packed stacking sequences of length N using group theory.

An algorithm has been developed that generates all of the nonequivalent closest-packed stacking sequences of length N. There are 2(N) + 2(-1)(N) different labels for closest-packed stacking sequences of length N using the standard A, B, C notation. These labels are generated using an ordered binary tree. As different labels can describe identical structures, we have derived a generalized symmetry group, Q approximately equal to D(N) x S(3), to sort these into crystallographic equivalence classes. This problem is shown to be a constrained version of the classic three-colored necklace problem.

Quantifying distortion from ideal closest-packing in a crystal structure with analysis and application.

A parameter, U(cp), that quantifies the distortion of the anion skeleton in a crystal from ideal closest-packing has been developed. It is a measure of the average isotropic displacement of the observed anions from their ideal equivalents. An ideal closest-packed structure can be fit to an observed structure by varying the radius of the ideal spheres, orientation and translation, such that U(cp) is minimized. Ideal structures were fit to the M1M2TO(4) polymorphs, pyroxenes and kyanite. The distortions of these crystals were analyzed in terms of the two parameters, U(cp) and the ideal radius. Changes in structures due to temperature, pressure and compositional effects were characterized in terms of these parameters.